STOCHASTIC CORRELATION MODEL 

OF GALACTIC BULGE VELOCITY 

DISPERSIONS AND CENTRAL BLACK 

g: HOLES MASSES 

o 

■ V. I. Dokuchaev? Yu. N. Eroshenko' 

^.^P | Institute for Nuclear Research of the Russian Academy of Sciences, Moscow 



> 
On 



X 



Abstract 



We consider the cosmological model in which a part of the Universe 
~ 10 -5 is in the form of primordial black holes with mass ~ 1O 5 M0. 
These primordial black holes would be centers for growing protogalax- 
ies which experienced multiple mergers with ordinary galaxies. This 
process of galaxies formation is accompanied by the merging of central 
black holes in the galactic nuclei. It is shown that recently discovered 
^ correlations between the central black holes and bulges of galaxies are 

5— ( ■ naturally reproduced in this scenario. 

CO 



1 Introduction 



Recent observations demonstrate that no less than 20% of regular galaxies 
contain the supermassive black holes (SMBHs) in their nuclei (Kormendy and 
Richstone, 1995; Ho, 1998). Different schemes were proposed for the origin 
of the central SMBHs due to: (i) gravitational instability and collapse of the 
supermassive star (e. g. Gurevich and Zybin, 1990; Lipunova, 1997); (ii) 
collapse of the dynamically evolving dense stellar cluster (e. g. Rees, 1984; 
Dokuchaev, 1991) or (iii) collapse of the central part of the massive gaseous disk 
(e. g. Eisenstein and Loeb, 1995). In all these scenarios the SMBHs are formed 
deep inside of the gravitational potential well of the galactic or protogalalctic 
nuclei. The quite specific possibility is the formation of primordial SMBHs in 
the early Universe (Zel'dovich and Novikov, 1967; Carr, 1975). 

There are observed three types of correlations between the mass of the 
central SMBH Mbh in the galactic nucleus and (i) stellar mass of galactic bulge 
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M b (Kormendy and Richstone, 1995; Ho, 1998), M BH ~ (0.003 - 0.006)M 6 ; 
(ii) bulge luminosity L B (Richstone, 1998), M B h/Lb — 1O~ 2 M /L , and (iii) 
velocity dispersion a e at the bulge half-optical-radius (Gebhardt et al., 2000): 

/ \ 3.75(±0.3) 

M M = 1.2(±0.2)X10»( 5S ^) M . (1) 

Correlations of types (i) and (ii) now are feebly marked and there is large data 
scattering with respect to the mean values of M B h- So recent observations 
provide only consistency between (i) and (ii) in spite of the absence of linear 
dependence between L B and M b . On the other hand the type (iii) correlation 
is more definite. 

The origin of discussed correlations is quite uncertain. The simplest as- 
sumption that the growth of the SMBH mass depends on the bulge processes 
meets with the problem of different scales: the galactic bulge scale is a few 
kpc, whereas the linear scale of accretion disk around of SMBH is much less 
than 1 pc. Some deterministic mechanism is needed for huge mass transfer 
from bulge to its innermost part. Silk and Rees (1998) proposed a feasible 
solution of this problem by considering the self-adjusting accretion flow from 
bulge under influence of radiation pressure on the early quasar phase of galac- 
tic evolution. However this model predicts the definite type of correlation, 
M B h oc <7g, with the power index differing from observed one. 

In this paper we explore the alternative approach by supposing that ob- 
served correlations are stochastic in origin. Our basic assumption is the ex- 
istence in the Universe of pre-galactic population of black holes with masses 
Mh ~ 1O 5 M near the recombination time. The similar hypothesis of the ex- 
istence of primordial massive black hole population was used by Fukugita and 
Turner (1996) for interpretation of quasar evolution. 

The supposed primordial BHs are mixed with dark matter due to there 
cosmological origin. So the total mass of these BHs in any galaxy Mh would 
be proportional to galactic dark matter halo mass M. As a result the corre- 
lation Mh oc M is primary in this model, but the aforementioned observed 
correlations M B h oc M b , L and cr" would be secondary and approximate in 
origin due to complicated process of galactic formation. We show below that 
requested relation M B h oc cr" with a power index a near to observed one fol- 
lows from our primary relation oc M with the value of a e determined 
by the dark matter halo mass M. The considered model reproduces also the 
correlation M BH oc M b but with a less accuracy. 

The other necessary requirement of our model is multiple merging of pri- 
mordial BHs with mass Mh into the one SMBH with mass M B h during the 
Hubble time. In Section [| we validate this requirement. It is known that for a 
single BH with mass -C 10 7 M Q the dynamical friction in the galactic halo 
is ineffective (Valtaoja and Valtonen 1989). Nevertheless for the early formed 
primordial BHs it is possible the process of dark matter "secondary accretion" . 
As a result the primordial BHs would be "enveloped" by the dark matter halo 
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with a mass of a typical dwarf galaxy and a steep density profile, p oc r ' . 
Indeed, the gravitationally bound objects formed at red-shifts z ~ 10 from the 
density fluctuations 5 ~ (at the recombination epoch). In the uniform 
Universe the BH with mass Mh ~ 1O 5 M produce this fluctuation inside the 
sphere containing the total mass Mh/S ~ 1O 8 M . We will call this combined 
spherical volume "BH +halo" by "induced halo" . The characteristics of these 
induced halo is defined in Sections [3] and |j. Our assumption of multiple merg- 
ing of primordial BHs may be violated in the galaxies of late Hubble types. In 
fact observations (Zalucci et al., 1998) demonstrate that masses of the central 
BHs in these galaxies are less than in E and SO ones and correlation (i) — (iii) 
are absent. 

Throughout this paper we consider for simplicity the flat cosmological 
model without the cosmological constant. We use index "i" for the quan- 
tities at the moment t j ~ 6 • 10 10 s of transition from radiation dominated 
to matter (dust) dominated phase. Respectively index "0" is used for quan- 
tities at the recent moment of time to- We call all BHs existing at nearby 
the recombination epoch by primordial ones. The other possible mechanisms 
of BH formation do not influence our results except the Section |^. We use 
term "bulge" as for elliptical galaxies and for central spheroidal parts of spiral 
galaxies which reminiscent the dwarf ellipticals. 



2 Primordial Black Holes 

Noncompact objects (NCOs) of mass ~ (0.1 4- 1)M & , consisting of weakly in- 
teracting nonbaryonic dark matter particles ~ (O.l-j-l)M like neutralino were 
proposed by Gurevich et al. (1997) for the explanation of microlensing events 
in Large Magellanic Clouds. The hypothesized NCOs or neutralino stars are 
originated from the cosmological fluctuations with a narrow sharp maximum 
~ 1 in the spectrum at some small scale. In addition to neutralino stars 
the same maximum in the spectrum of cosmological fluctuations produce also 
the massive primordial BHs with mass ~ 10 5 M Q (Dokuchaev and Eroshenko, 
2001). So the hypothesized dark matter NCO and primordial BHs may be 
indirectly connected through their common origin from the same cosmologi- 
cal fluctuations. The spectrum with a sharp maximum at some scale arises 
in some inflation models (Starobinsky, 1992; Ivanov et al., 1994; Yokoyama, 
1995; Garcia-Bellido, 1996). At the same time the spectrum beyond the maxi- 
mum may be of the standard Harrison-Zel'dovich form and reproduce the usual 
scenario of large-scale structure formation in the galactic distribution. 

Adiabatic density fluctuations of the matter on the scale less than the 
horizon grow logarithmicallly in time during the radiation domination epoch. 
These fluctuations with amplitude 8i ~ 1 at the moment U corresponds to 
radiation density fluctuations ~ 0.05 at the moment th when disturbed region 
enter into the horizon (Dokuchaev and Eroshenko, 2001). So large fluctuations 
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at radiation dominated epoch produce BHs near the time th (Zel'dovich and 
Novikov, 1967; Carr, 1975). 

Formation of BHs takes place on the tail of the Gaussian fluctuation distri- 
bution, whereas the main part of NCO produced from the r.m.s. fluctuations. 
Therefore only small part of fluctuations which result in formation of NCO 
may produce BHs at the radiation domination epoch. In other words due 
to the large threshold of BH formation the major part of fluctuations do not 
collapse into the BHs and evolve continuously up to the end of radiation dom- 
ination epoch. At the radiation domination epoch the mass in the comoving 
volume varies as M(t) = M x a(U) / a(t) , where the scale- factor of the Universe 
a oc t 1 / 2 and M x is the comoving mass at the moment of transition to the 
matter domination. The mass M x equals approximately to the mass of non- 
relativistic matter inside the fluctuation, i. e. the mass of NCO which may be 
formed from this fluctuation. On the other hand 

Att 

M(t) = -(2ctfp(t), (2) 

where p(t) = 3/32irGt 2 . From these relations for M(t) we find the time t h and 
mass M h of BH formation in dependence of the mass of NCO: 

tk = ( M, G ,^3 =1 ,g) 2/3 (_|_) 1 ' 3 3, (3) 

M h = cMfG-"f J = 2.3 ■ 10» ( jr A ;r -) 1/3 M e . (4) 

In our model the NCOs add BHs are originated from the fluctuations of the 
same type but at different moment of time. Whereas the strong difference 
in masses of NCO add BH is accounted by the mass variation in time inside 
the fixed comoving volume: the mass of radiation at the radiation domination 
epoch (e. g. at moment t ~ 1 s) far exceeds the cold dark matter mass at the 
matter domination epoch in the same comoving volume. 

There are definite astrophysical limitations on the number and mass of 
primordial BHs: 

1. From the Universe age limit it follows that fraction of BHs < 1. Addi- 
tionally the PBHs with Qh ~ 1 would disturb the background microwave 
spectrum if they formed after the time ~ 1 s of electron-positron pair 
annihilations (Carr, 1975). 

2. The possibility of tidal destruction of globular clusters by the primordial 
BHs gives the BH mass limit M h < 1O 4 M if these BHs provide the 
major contribution to the dark matter (Moore, 1993). 

3. The model-dependent calculations of the contribution from matter accre- 
tion by PBHs at pre-galactic and recent epoch into background radiation 
give approximately tt h < 1(T 3 -=- KT 1 for M h ~ 1O 5 M (Carr, 1979). 
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4. From the absence of reliable events of gamma-ray bursts (GRBs) lensing 
it is obtained the limit Qh < 0.1 for PBHs in the intergalactic medium 
with masses 1O 5 M < M h < 1O 9 M (Nemiroff et al., 2001). An even 
more stringent limit tt h < 0.01 for 1O 6 M < M h < 10 8 M Q is obtained 
from VLBI observations of the lensing of compact radio-sources (Wilkin- 
son et al, 2001). 

In this paper we consider the case Qh ~ 10~ 5 in accordance with all aforemen- 
tioned limits. 



3 Formation of Induced Halo 

A massive induced halo (IH) or a heavy dark matter envelope around the PBH 
is formed due to dark matter accretion in the expanding Universe. Spherical 
matter accretion onto the compact object (SMBH or galaxy) produces the 
stationary matter density profile of the form p oc r~ 9//4 (e. g. Ryan, 1972; 
Gunn, 1977). Here we reproduce the similar calculations in the form suitable 
for binding the mass of PBH with a mass of gravitationally trapped matter 
and a corresponding red-shift. 

The matter in the spherical layer of radius r around the PBH in the uniform 
Universe exerts the acceleration 

d 2 r G /lir Air , . 

where and pi = 4.4 • 10~ 18 (/i/0.6) 8 g cm -3 is correspondingly the radius of 
the spherical layer and the Universe density at the moment £, when the strong 
growth of fluctuations is beginning, and h is the Hubble constant in units of 
100 km s" 1 Mpc -1 . Let us multiply (|j) on db/dt and integrate it by using 
normalization r(t) = Tib(ri,t) and initial condition b(U) = 1, b(ti) = H(ti)b. 
After using an additional normalization r = (t — t^a 1 ^ 2 , where 

(2GM h 8ttG \ 

the integrated (H) takes the form: 

b{db/dr) 2 = l-bE, (7) 

with E = (a — H 2 (ti))/a. Equation of this type admit an exact solution 
(Saslaw, 1989). From this solution we find the expansion termination moment 
of the considered layer: 

Ufa— U. 8 



* 4 \M h 



Here M s = 47rpjr 3 /3 is the mass inside the layer with the exception of central 
PBH mass. At the moment of expansion termination b(t s ) = 1/E and sub- 
sequently the layer is contracted up to the moment t co i ~ 2t s . From (ffl) by 
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using relation t = £o/(l + z) 3 ^ 2 for the flat Universe without the cosmological 
constant we obtain 

z col « 0.36^—^ - 1 « 2 • 10 3 — p - 1, 9 

where z< = 2.4 • 10 4 /i 2 . 

We will suppose that spherical layer is detached from cosmological ex- 
pansion and virialized and after its expansion termination and a subsequent 
contraction from radius 7*j&(rj, t s ) to radius r co i = rib{r i: t s ) /2. For this virial- 
ization radius of the IH around the central PBH we find 

r- (J-)'" = 0.45 L^f ( » do) 



2£ 2M h ^4^; V 105M o/ Vl + ^oi 

Now we obtain the requested relation for the mass and radius of IH 

M, ( r^ = ( 2M ^)^(^) 1/4 

= 3.7- 10' (^-) m (^L-) m (±X M B (U) 

lio 5 M y v lk P c / V - 6 / 



and density distribution in the IH 

1 dM s (r) 



p(r) 



Attt 2 dr 



3 / 4 / „ \ - 9 / 4 / u \ 2 



1.7.10-{-^. ^- A] gcm -3. (12) 

llO 5 M i V lk P c / V - 6 / 



4 Protogalaxies Growth Termination 

The growth of IH terminates at the epoch of nonlinear growth of ambient den- 
sity fluctuations with mass M and radius R producing the same gravitational 
acceleration: 

GM s (r col ) = GM 
(2r col ) 2 (2R)*' [ } 

where radiuses 2r co ; and 2R correspond to the moment of expansion termina- 
tion. At this condition the capture of new spherical layers is shutting down 
(see. (H)). 

The r.m.s. fluctuation in the mass scale M is defined by expression (Bardeen 
et al, 1986) 

° 2 (M, z) = * + J dkk 2 P (k)W 2 (k } M). (14) 

Here W(k, M) is a filtering function and P(k) is a power spectrum: 

(S*M = (2n) 3 P(k)S ( §\k - k'), 5 h = j S(r)e lkr d\ (15) 
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where 8$ (k — k') is the Dirac delta-function and angular brackets denote the 
ensemble averaging. Let us choose Po(k) for the model with a cold dark matter 
(Bardeen et al., 1986) with the Harrison-Zel'dovich initial spectrum: 

Po(k) = aHn ( 7 3 4^ 2 34 ^ (16) 

where q = kVL^ l h~ 2 is written in the comoving coordinates with units Mpc -1 . 
For simplicity we put the dark matter density parameter = 1. Normaliza- 
tion constant a is chosen from the condition oo = 1 at the scale 8/i _1 Mpc. 
The condition for object formation with the mass M at the red-shift z co \ is 

a(M,z col ) = 8 c , (17) 

where for the model of spherical contraction 8 C = 3(37r/2) 2/,3 /5 (White, 1994). 
Radius of ordinary protogalaxy (without the PBH) after its cosmological ex- 
pansion termination and virialization is (Saslaw, 1989; White, 1994) 

D 3R t (3M\ 1/3 
/i — - - - , Jx 



105i' \4n Pl J y ' 

By taking into account ([K]), flTFD , (Jl8|), and putting Si = &i(M) the condition 
(P~3|) can be rewritten in th form: 



\ 2 J \5j l + z col ) 



:,,,/) ^ 

It follows that M s = M, because in ( |10D and (^8|) it is supposed that IHs and 
ordinary protogalaxies have the same coefficient of nonlinear contraction equals 
to 0.5. We solve numerically the system of equations fllTf) and ( |I9| ) relative to 
independent variables z co \ and M with the M h as parameter. See Fig. 1 for 
calculated relations for z col (M h ) and M(M h ) = M s (M h ). For M h = 2.3-lO 5 M 
we find z co i = 8.8 and M = 7.2 • 10 7 M Q . As a result up to epoch z = 8.8 the 
PBHs with mass Mh = 2.3 • 10 5 M Q had time to capture an additional mass 
which ~ 300 times exceeds the PBH mass. 



5 Merging of Black Holes 

In the preceding Chapter we demonstrate that massive IHs with mass M s = 
7.2 • 1O 7 M are formed around the PBHs. These IHs are massive enough to 
sink down to the galactic center during the Hubble time under influence of 
dynamical friction. Though the fate of nested PBHs inside the central parsec 
of the host galaxy is rather uncertain. Valtaoja and Valtonen (1989) considered 
interaction of central BHs after merging of galaxies. The late phase of two BHs 
merging in the galactic center depends on many factors (Menou et al., 2001). 
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Without the detailed elaboration we will follow to Menou et al. ( 2001) by 
supposing that multiple PBHs merge into a single SMBH during the Hubble 
time. 

Notice that density (|i~2|) strongly grows towards the center and smoothed 



out only at the distance where M s = Mh according (11). For Mh = 
2.3 • 10 5 M Q we obtain ~ 1 pc and density p ~ 1O 4 M ~ 3 . By using Chan- 
drasekhar time for dynamical friction (see e. g. Saslaw, 1989) of the PBH with 
mass Mh = 2.3 • 1O 5 M we obtain the estimation of characteristic time for 
PBH merging 

( '~4^W~ 5 ' 105yrS ' < 2 °> 

where v ~ (GM h /r h ) 1/2 is the PBH velocity, A « 10, B « 0.426. The PBH 
of larger mass would be merged even faster. As a result the late phase of 
PBHs merging lasts very fast, and the probability simultaneous presence in 
the galactic nucleus of three or more BHs is low. 

Coalescence of PBHs in the galaxies must be accompanied by the strong 
burst of gravitational radiation. The projected interferometric detector LISA 
is capable to detect the coalescence events up to the red-shiftsz ~ 10, if BH 
masses no less than 1O 3 M (Menou et al., 2001). For the calculation of grav- 
itational bursts distribution there is needed the numerical simulation of hi- 
erarchical BH and galactic merging. Menou et al., (2001) carried out the 
corresponding simulation for ordinary galaxies. A simple estimation of the 
burst rate from the observable Universe gives 

where n g is a mean number density of structured galaxies and AMs a mean 
number of merging per galaxy. So there is a principal possibility for the veri- 
fication of considered model by the LISA detector. 



6 Origin of Correlations 



Fluctuation spectrum (jig) in the confined mass region can be approximated 
by the power law with the effective index n = dlnPo(k)/dlnk. According 
to the formation condition fll"T|) for power law spectrum the effective galactic 
mass formed at red-shift z is (White, 1994) 

M = M (l+zy^3, (22) 

where M = const. For galactic mass M = 1O 1O M we obtain n = —2.28 
M = 2.5 ■ 1O 16 M and for M = 1O 12 M we find correspondingly n = -1.98 
M = 7 ■ 1O 14 M . Velocity dispersion is estimated as 

9 GM ( 3M \ 1/3 . . , . o 
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For PBH cosmological density parameter Qh and effective PBH merging in the 
galaxies the preceding relation gives the final mass of the central BH 



_L2_ _n + 3 /4 7r G , 3 po \ l-r 



M BH =ipn h M= ipn h at n M {—^J > ( 24 ) 

where factor ijj is responsible for the possible additional growth of the central 
BH due to accretion. From (1241) with h = 0.6 it folows 



(3.66^-4.03) 

M BH = (1.1 + 1.7) -ltf I ( — p _J M , (25) 




where pair of coefficients 1.1 and 3.66 corresponds to M — 10 10 M Q , and re- 
spectively pair 1.7 and 4.03 correspond to M = 1O 12 M . In summary the 
considered model is in a good agreement with observation data ([!]). The fluc- 
tuation spectrum at the galactic scale, n « —2, completely defines the power 
index a ~ 4 in the relation Mbh a®. It is easily verified that for Qh ~ 10 -5 
the contribution IHs (there mass is determined in Section |]) to the total galac- 
tic mass is negligible. 

Correlation of the form Mbh oc Mb is explained by the common relation of 
dark matter mass in galaxies of all types with their total baryonic (star) mass 
with the main input from old stars. In the framework of hierarchical model 
of galaxy formation the recent galaxies formed from multiple merging of low 
mass protogalaxies with early star formation proceeding, in this case the old 
star population of spherical subsystem is protogalactic in origin. 

The part of galactic mass contained in stars becomes a fixed and indepen- 
dent part of the total galactic mass M due to statistical averaging after the 
large number of protogalactic merging. So the total star mass in the galaxy is 
Mf, = fj,f s M, where m 0.05 is the baryonic mass fraction of the Universe, f s 
is the part of baryons passing into the stars. After sinking down to the galactic 
center and merging of all PBH the resulting mass of the BH is 

(26) 

which with the order of magnitude corresponds to observations Mbh — (0.003— 
0.006)M b . 

The central BH masses in Sa, Sb, Sc galaxies are less in a mean than in 
E and SO galaxies (Salucci et al., 1998). In our model this is related with a 
relatively late formation of Sa, Sb, Sc galaxies when the main part of PBHs 
do not have enough time to sinking to the galactic center. In particular ~ 10 2 
PBHs of mass M^ ~ 1O 5 M can inhabit our Galaxy. 
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7 Conclusion 



We explain the observed correlations between bulge parameters and central 
BH masses in galaxies by the multiple merging of PBHs of mass ~ 1O 5 M 
generated in the early Universe. This model predicts the existence of BHs 
with mass ~ 10 5 M Q in beyond the dynamical centers of spiral galaxies and 
in the intergalactic medium. Possibly one of these BHs was detected by the 
Chandra Observatory in the galaxy M82 (Kaaret et al., 2000). The observa- 
tional signature for verification of this model is the gravitational burst after 
BH merging in the galactic nuclei which may be detected by the projected 
laser interferometric detectors of the LISA type. 
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Russian Foundation for Basic Research grants 01-02-17829, 00-15-96697 and 
00-15-96632. We also acknowledge the referees of this paper for helpful re- 
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FIGURE CAPTIONS 

Fig. 1. Functions z co i(Mh) and M(Mh) from the numerical calculation of 
the system of equations (|PT|) and (|19|). 
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